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Small regulatory RNA (sRNA) that acts by an antisense mechanism is 
critical for gene regulation at the posttranscriptional level. Recently, an Hfq- 
dependent sRNA named MicM, which is related to the regulation of outer 
membrane protein, was verified as a novel antisense sRNA due to its 
catalytic mode of regulation. Here we propose a simple kinetic model for 
the enzyme-like regulation mode of sRNA and study in detail the noise 
properties of the target gene under various recycling rates of the regulator. 
Keywords: We predict that the recycling rate of sRNA and other relative parameters 
small RNA; have significant influence on the noise strength of target expression. In 
noise; comparison with the stoichiometric regulatory mode, a lesser fluctuation of 
RNA trap target expression was observed near the threshold at which the transcrip- 
tion rates of both sRNA and target mRNA equal each other. We also found 
that the new mode is better in terms of rapid response to external signals. 
However, it needs more time to achieve target recovery if the stimulating 
signal disappears. Additionally, the obtained time evolution results of the 
MicM-ybfM interaction system based on our model are consistent with 
previous experimental results, serving as experimental evidence to back up 
our theoretical analysis. 

© 2010 Elsevier Ltd. All rights reserved. 


Numerous studies sought to understand the 
dynamic features of gene expression controlled by 


Introduction 


Small regulatory RNA (sRNA) is a type of 
noncoding RNA with typical length varying from 
~50 nucleotides to several hundred nucleotides. As 
a crucial regulator, it is involved in many biological 
functions such as stress response, * quorum 
sensing,” ” and so on. Approximately 80 sRNAs 
have been verified® in Escherichia coli, most of which 
use an antisense mechanism by base-pairing with 
other RNAs.' These sRNAs are often induced under 
stress conditions and respond rapidly to environ- 
mental threat.” '° They are usually retained at a very 
low copy number, so stochastic effects can become 
prominent. 


*Corresponding author. E-mail address: shihl@itp.ac.cn. 
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Abbreviations used: sRNA, small regulatory RNA; 
SSA, stochastic simulation algorithm; MFPT, mean first 
passage time. 


sRNAs that lead to coupled degradation.!!-'4 
Comparison of protein and sRNA regulators based 
on simulations drew the conclusion that regulation 
by sRNA is advantageous when a rapid response to 
an external signal is required. This is in consistent 
with experimental data.!*+ Noise analysis indicat- 
ed that the noncatalytic mode of gene regulation by 
sRNA has a unique feature in noise resistance. 12-14 
Different from the ordinary mechanism of base- 
pairing sRNA regulators that act stoichiometrically 
on their targets,! recent studies found that, in 
Salmonella, an SRNA named MicM leads to selective 
degradation of the sRNA-mRNA duplex.!>!” Its 
binding to chb mRNA leads to coupled degradation 
of both sRNA and mRNA. For another target-ybfM 
mRNA, the sRNA part of the RNA duplex remained 
undestroyed. However, the mRNA part could not 
escape the fate of degradation. Experimental results 
indicated that MicM-mediated regulation requires 
Hfq protein to enhance RNA duplex formation”? 


0022-2836/$ - see front matter © 2010 Elsevier Ltd. All rights reserved. 
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and implied that one molecule of MicM can 
participate in the degradation of several molecules 
of ybfM mRNA.” 

The catalytic characteristic of SRNA enriched our 
understanding of gene regulation: It is a fusion of the 
traditional mechanism of the sRNA regulator and 
enzyme catalysis. It is natural to ask whether this 
new mode of posttranscriptional gene regulation 
mediated by sRNAs is well suited for biological tasks 
and whether this mode of gene silencing is efficient 
in suppressing fluctuations in its target genes. 

Here we propose a simple kinetic model for 
catalytic-sRNA-mediated gene silencing to investi- 
gate the noise properties of a target gene and the 
time evolution of the interaction system by Monte 
Carlo simulations. Some unusual features of this 
regulatory mode are found and compared with 
ordinary antisense mechanism. Catalytic sRNA is 
promising as a more effective regulator owing to its 
fast response ability and low recovery capability of 
target expression. 


Results 


A basic model for catalytic-sRNA-mediated gene 
silencing 


Our theoretical study of catalytic-sRNA-mediat- 
ed repression starts with a mass action model 
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(Fig. 1). The minimal system consists of a single 
species of small RNA (s), a single species of 
mRNA target (m), the duplex (s_m) formed by the 
two species of RNA, and the protein product (p) 
of the mRNA. As a simplified model for catalytic 
sRNA, an sRNA can bind the 5’ untranslated 
region of an mRNA at its ribosome binding site to 
block translation initiation. In this model, we 
assume the following: 


(1) Small RNA and its target mRNA are coded by 
independent genes using different promoters for 
their transcription. They are synthesized at rates 
a, and Am, respectively. 

(2) The binding of sRNA to mRNA is irreversible 
and fast with the binding rate k, while an sRNA- 
mRNA duplex is degraded with the linear 
rate Bsm. A proportion of sRNA is successfully 
recycled as provider of sRNA regulators from 
the degraded duplex; the recycling ratio ô ranges 
from 0 to 1. 

(3) Proteins are translated at a linear rate k, off 
mRNA molecules, and each mRNA molecule 
produces on average of n protein molecules 
during its lifetime = so k,="Bm. 

(4) All molecules of sRNAs, mRNAs, sRNA- 
mRNA duplexes, and proteins, are subject to 
degradation with linear rates Bs Bm Bsm and Bp, 
respectively. 
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Fig. 1. Schematic diagram of catalytic-sRNA-mediated gene regulation. 
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Steady state solutions(nM) 
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The process, as shown schematically in Fig. 1, 
can be quantitatively described with a system of 
differential equations: 


ds 
dt = as — Bss — ksm + öBsm(s-m) 
dm 
= Xm — Bum — ksm 
dt (1) 
d(s-m) 
dt, = ksm — Psm (s-m) = Bga(s-m) 
P 
di = kpm = Byp 


Here, we assumed that the other proteins involved 
(Hfq, degradosome, etc.) are present at levels that 
are sufficient not to limit reaction rates. The 
degradation of all molecules is caused by direct 
degradation and growth dilution. Direct degrada- 
tion of the sRNA-mRNA duplex will save a 
fraction of sRNAs, but growth dilution will not. 
In order to show this difference, we explicitly write 
the growth dilution term for the sRNA-mRNA 
duplex, where Bga is related to the growth rate. 
Usually, direct degradation of the sRNA-mRNA 
duplex is much faster than growth dilution, and 
we neglect it in the ensuing discussion. The 
steady-state solutions for s and m can be deduced 
from the following equations: 


a; — Bss — k(1— 5)sm = 0 
Am — By, — ksm = 0 


(2) 


For mRNA, the rate k here is the same as the 
coupled degradation rate in the previous stoichio- 
metric mode” (Fig. 2a); however, for sRNA, part 
of the sRNAs (dksm) can be recycled from the 


Fig. 2. Steady states with different recycling ratios. (a—c) 
Steady-state solutions of mRNA (left ordinate), sRNA, and 
sRNA-mRNA (right ordinate) with 5=0, 0.5, and 1, 
respectively. (d) The levels of target mRNA in steady 
state with different recycling ratios of sRNA. The lines are 
obtained from steady-state solutions (Eqs. (3) and (4)). The 
symbols are obtained from Monte Carlo simulations, 
where 10° steps of simulation were used for each recycling 
ratio 5. Parameters are set as a,=1 nM min ~~ and 
k=0.1 (nM min) '. The other parameters are estimated 
from experiments.” We set B.=By,=0.02 min™! and 
Bsm=0.1 (nM min)~’. The half-lives of both sRNA and 
mRNA are about 27 min, and the half-life of the RNA 
duplex is about 4 min." 


198 


Analysis of Catalytic- sRNA-Mediated Gene Silencing 


degraded RNA duplex. Thus, the solutions under 
reaction equilibrium conditions are (for 6€[0,1)): 


Fe | EE. 
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* = 2 m* 3 
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where N= is the leakage rate of mRNA 
expression”? (Fig. 2a). The effect of the leakage 
expression rate of mRNA can be neglected, since 
the binding reaction is fast due to the participation 
of Hfq. In this study, k is set as 0.1 (nM min)! to 
represent a fast combination between sRNA and 
mRNA. The catalytic mechanism maintains the 
characteristic of a “threshold-linear response.” 1°70 
The threshold here, however, is characterized by 
74%, not by a,. It rises as 5 increases. Above the 
threshold, the mRNA level will start to increase 
linearly with this production rate a. It is obvious 
that at a fixed value of as, a higher recycling ratio of 
sRNA leads to a higher level of gene silencing (Fig. 2). 

The sRNA of the degraded RNA duplex will be 
100% recycled if 6=1, and the recycled sRNA will 
exert its effect again on other mRNA molecules. This 
process is similar to a protein regulator taking part 
in posttranscriptional regulation, which results in 
the degradation of target mRNAs. Taking limit 6— 1 
in Eq. (3) or deriving directly, we have: 


st = at 
Bs 
m* = a = ae (4) 
1+ La 1+ ae 
Bin ` 


where mọ = “ denotes the stationary level of 
mRNA in the 'absence of the regulator. Thus, the 
level of sRNA molecules remains constant even 
when Qn is varied, whereas the number of sRNA- 
mRNA duplex linearly increases with a,, (Fig. 2c). 
With parameters estimated based on experiments, 
mRNA is suppressed at a low level even when 
Qn >> as (Fig. 2d). It is interesting that the steady 
state of this mode is similar to that of protein-based 
transcriptional inhibition, which is characterized by 


m* = T E" Here r is the concentration of the 
a 


transcription factor, a is related to binding affinity, 
and n is the Hill coefficient.™™?? It should be 
emphasized that, in both cases, the efficiency of 
regulators is closely related not only to the 
concentration of the regulator but also to the 
binding rate with its target. Catalytic sRNA inter- 
polates between the conventional protein regulator 
and the sRNA regulator. 


Description of noise strength 


It is well known that gene regulation is intrinsi- 
cally noisy due to the stochastic nature of underly- 
ing biochemical reactions. °- In this article, we 
define intrinsic noise as fluctuations in mRNA 
number and output protein number at a fixed 
input of Qm, Qs, and so on. 

The master equation for the distribution of the 
mRNA, sRNA, sRNA-mRNA duplex, and protein 
copies (m, s, c, and p, respectively) in this model is 
given by: 


Š Pm, 5,7.) = fen (ha a 1) E Br (h T 1)m 


+ &s(h - 1) + B, (hf - 1)s 
+ k(hihth. - 1)ms 


m s € 


+ Bsn (A, he ~ 1)c + Psm (he ~ 1)c 


+ kpm(h, = 1) 
+ Bp (hy = 1)p P(m,s,p,c) (5) 


where h* and h~ are, respectively, the lowering and 
raising operators defined through hyf(x)=f(x+1), 
hyf(x)=fx—-1), and hxf(0)=0 for any function f(x). 
All time derivatives are set to zero here, as we focus 
on the steady-state properties of gene expression. 

To describe noise strength, we use a sensitive 
measure of noise, the Fano factor, which is defined 
as Fy = = for stochastic variable x, where o2 is the 
variance and yu, is the mean. For the Poisson process, 
the factor equals 1. 


Dependence of noise strength on recycling ratio 


It has been previously found that sRNA-based 
regulation is more efficient than transcription-factor- 
based regulation both in responding rapidly to 
signals and in filtering out noise from the 
environment.'**° Here we focus on the recycling 
rate of the sRNA regulator in the catalytic mecha- 
nism and quantitatively describe its influence on the 
noise strength of target expression. 

It is difficult to directly solve the master equation 
describing the intrinsic stochastic behavior of small- 
RNA-mediated regulation, so we use the well-known 
stochastic simulation algorithm (SSA) Gillespie direct 
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method*® to generate the stochastic trajectory of 
the regulation process and then analyze the 
characteristic of noise (see details in Materials 
and Methods). 

The threshold-linear mode of gene expression 
regulated by ordinary antisense sRNA has been 
carefully analyzed. It has been concluded that the 
fluctuation of target gene expression depends 
largely on the ratio of the transcription rate of the 
target mRNA (am) to the transcription rate of the 
sRNA regulator (as). 77 A similar behavior is 
now also predicted for gene expression controlled 
by catalytic sR NA. Based on the master equation 
(Eq. (5)), the average number of mRNA or the first 
moment of m was drawn in Fig. 2d. The result shows 
that the average numbers of species provided by 
stochastic simulations are consistent with the solu- 
tions of the deterministic equation. 

We present the relationship between the Fano 
factors of mRNA and the production rate a,,, under 
different recycling ratios in Fig. 3a. It shows that the 
inflection points of these curves move rightward 
along the a,, axis. The peak position changes as 
y} %s, and the peak value becomes larger as ô 
increases. Note that in the regime where the 
transcription rates of mRNA and sRNA are compa- 
rable, the value of the Fano factor decreases as 6 
increases. This result may be important in biological 
function, as is mentioned in the literature (that most 
actual biological parameters are near this point 
when sRNA silences its target gene). In other words, 
the catalytic-sRNA-mediated regulation mechanism 
tends to attenuate noise more than ordinary sRNA 
regulation when gene silencing happens. Due to the 
large stochastic fluctuation in the cell, it is important 
to ensure that the silencing function is reliable. 

Protein noise has a similar behavior, except for 
%m>>z+s% (Fig. 3b), where the Fano factor of 
protein approaches 1+ - 2 "y 


(equal to about 2 


Bin 
here) instead of 1. The same result can be deduced 
from the master equation in the following manner: 
Multiplying the master equation by m, m?, s, and 32, 
respectively, and summing over m, s, p, and c, with 
some approximations, we can derive that for 


Am K 7h and dm> h %s, Fm indeed approaches 
1. It was found that Fpl for Am THis and that 


= b 1 _ kh ; 
Fp=~1+ ne for Am 75%, Where b = 3, is the 


average number of proteins produced from one 
transcript (called “burst size” by some authors).*”*? 
The “burst” of protein in an mRNA translation 
process causes amplification of noise in gene 
expression, which has been elucidated in detail.7°” 

A previous work on stoichiometric sRNA-medi- 
ated gene regulation” found that there is bimodal- 
ity of target mRNA under certain parameters near 
the threshold, and the high correlation between 
sRNA and its target mRNA will form an effective 


Om (NM/min) 


Fig. 3. Fano factors for target mRNA (a) and protein 
(b) with different recycling ratios. The gray broken line 
in (a) represents the Poissonian limit F,,=1, and that in 
(b) represents the limit 1 + a Bp is set to 0.02, and the 


T 
other parameters are set as in'Fig. 2. These results are 
obtained from 10° steps of Monte Carlo simulations. 


potential barrier to shape the bimodal distribution. 
We also found bimodality in this catalytic-sRNA- 
mediated regulation mechanism, which needs to be 
tested by future experiments (Fig. 4). 


Comparison between stoichiometric mode and 
full catalytic mode 


We then focus on two limit cases in our model 
and find out their different behaviors. In our model 
(Eq. (1)), 5=0 and 6=1 correspond to a stoichio- 
metric mode” and a fully catalytic mode, respec- 
tively. In a real biological system, degradation of 
the sRNA-mRNA duplex is caused by direct 
degradation and growth dilution. Direct degrada- 
tion will save a fraction of sRNAs, but growth 
dilution will not. Thus, the recycling ratio of sRNAs 
will never exactly equal 1. But direct degradation is 
usually much faster than growth dilution, and 6 is 
possibly close to 1. Thus, it is still helpful to study 
the behavior of the system at the limit case of 5=1. 
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Fig. 4. Bimodal distribution of target mRNA. The 
parameters are set as: a,=1 nM min +, Q,=3 nM min /, 
Bs=Bm=0.02 min™!, B,,,=0.1 (nM min) !, X\=0.0001 nM 
min |, and 6=0.6. The ordinate is scaled by the total 
number of mRNA molecules. 


As 5=0, it becomes the previous model that was 
studied in detail.” When a,,<a,, almost all 
mRNA molecules will be targeted by small RNA, 
and only a thimbleful of mRNA, whose amount is 
not correlative to small RNA fluctuation, leak out. 
Alternately, when a„ > as, small RNA is used up, 
and the majority of mRNA is alive for translation. 
These two situations are close to the Poisson 
process, which leads the Fano factors' approach 
to 1 (Fig. 5, black line). 

In the case of 6=1, our result shows that the Fano 
factor remains close to 1, although the mRNA 
production rate (a,,) changes (Fig. 5). When 
Qn <Q, the result is similar to Levine's model,” 
where mRNA molecules are mostly degraded, 
resulting in its fluctuation being irrelevant to small 
RNA fluctuation. While a,,>a;, over quite a large 
range, mRNA will be consumed, as sRNA can 
remain at a constant level. The sRNA here can be 
recycled like an enzyme after the degradation of the 
mRNA-sRNA duplex. According to our simulations, 
Fm remains close to below 1 even when am > as. 

To obtain a better understanding of full catalytic 
mode, we will return to deterministic results (Eq. (4) 
and Fig. 2c). It is obvious that more and more 
mRNA molecules will bind to sRNA molecules and 
then degrade as a,n increases. However, the level of 
free sRNA molecules stays the same because the 
effects of binding to mRNA and recycling from the 
duplex will cancel out each other. Small RNA 
pairing with its target mRNA produces a double 
effect: enhanced degradation of mRNA and protec- 
tion from direct degradation (Eq. (4)). The latter 
effect enables sRNA to always satisfy the need for 
silencing by adjusting the level of the mRNA-sRNA 
duplex. That is, when 5=1, although free sRNAs 
stay at the same level, the total SRNA concentration 
does increase with am, and this is why mRNAs are 


always repressed by sRNA and mRNA noise 
remains at a quite low level. 


Influence of other parameters on noise 


We analyzed the influence of the sRNA recycling 
rate on noise and found that a strong fluctuation of 
gene expression can be produced by a large 6. With 
fixed transcription rates (a, and a,,), aS more 
regulators are recycled, more targets will be 
influenced. It must be noted that sRNA recycling 
only happens with the degradation of the sRNA- 
mRNA duplex. Thus, the rate parameters of SRNA- 
mRNA duplex formation (k) and degradation (Psm) 
should be considered as well. Here we give a 
quantitative view of how these parameters affect the 
noise of both mRNA and its protein product. 

Figure 6a shows that noise is enhanced when the 
degradation of the duplex is enhanced. For certain 
recycling ratios, a fast degradation will speed SRNA 
recycling without changing the inflexion of the Fano 
factor curve, which is only decided by the recycling 
ratio. An upper limit of Bsm, above which sRNA 
recycling rates will not increase further, is expected. 
When PBsm=0, the model reduces to an ordinary 
sRNA-mediated model.” 

Figure 6b shows how the binding rate k influences 
the fluctuation of output protein. We calculate the 
Fano factor of protein at a different binding rate k for 
5=0.5. Over a range of &m< as, the Fano factor of the 
output protein decreases with increasing k. There 
are two limits of F, as the binding rate (k) is varied 
over a large range (from 0.0001 to 10 (nM min)~’). 
When the process of binding is very fast, the few 
proteins produced for mRNA are fully repressed, 
and the noise curve approaches the Poissonian limit. 
For k~0, mRNA is freely translated, and the Fano 
factor approaches the limit of 1+b (=2). In the 
silencing region, sRNA-based regulation can filter 
noise in the mRNA.’ A larger binding rate k means 
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Fig. 5. Fano factors of mRNA for stoichiometric mode 
(black continuous line) and catalytic mode (gray continuous 
line). The light-gray broken line represents the Poissonian 
limit Fm=1. 
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Fig. 6. Influence of Psm and k on noise at a fixed 
recycling ratio 5=0.5. (a) Fm is affected by the degradation 
rate Bsm of the sRNA-mRNA duplex. (b) F, is affected by 
the binding rate k. Here a,=1 nM min™!. All other 
parameters are set as in Fig. 2. 


a more efficient sRNA, which leads to a smaller 
fluctuation of mRNA in the silencing region. 

On the contrary, the result shows that over a range 
of Qm” Qs, a larger binding rate k of mRNA with 
sRNA produces a larger fluctuation of the output 
protein. At this region of a,,,, the noise of sRNA will 
propagate to the mRNA, and the efficiency of noise 
propagation depends on the binding rate. As k 
approaches infinity, the number of free mRNA that 
can work for translation is the difference between 
the number of mRNA and the number of sRNA. 
Thus, it will keep all noise originating from mRNA 
and sRNA in the free mRNA. With a decreasing 
combination rate k, the coupling between mRNA 
and sRNA becomes weaker, and the noise of mRNA 
propagated from sRNA becomes smaller. 


Response and recovery of target expression 
It is important for a cell to be sensitive and to 


rapidly respond to an external signal. For the new 
sRNA regulation mechanism, the response time 


characterizes how long it takes for mRNA to change 
from an sRNA-free state to an sRNA-bound state. 
Below, we will use the mean first passage time 
(MFPT) to measure this important characteristic time. 

Because of stochastic reactions and molecular 
fluctuation, it is necessary to calculate the time it 
takes to reach a new state in probabilistic meaning. 
The previous work provided information on re- 
sponse time through experiments.'? However, 
solving the master equation to obtain an analytic 
result for response time is difficult. Here we still use 
the Gillespie direct method to solve the problem 
concerning both stochastic reaction and molecular 
fluctuation. From the results of Fig. 7a, we see that a 
faster response of the target always corresponds to a 
higher recycling ratio of the regulator. A smaller 
response time means a more sensitive and effective 
regulatory mode. 

Recovery time is the time required for output 
recovery after the external signal has been removed. 
It is calculated similarly as the response time (see 
Fig. 7b). As the recycling ratio increases, the target 
mRNA will decrease, and more time will be 
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Fig. 7. Response time (a) and recovery time (b) of target 
gene expression. In the simulations, the threshold is set to 
70% of the difference between the initial-state mRNA level 
and the final-state mRNA level. a, is set to 0.5 nM min “1, 
and the other parameters are set as in Fig. 2. 


202 


Analysis of Catalytic- sRNA-Mediated Gene Silencing 


required for the recovery of this target gene 
expression. In the limit case of 5=1, the level of 
protein cannot recover from the repressed state even 
in several cell cycles (Fig. 7b). However, if the 
binding of the two RNA molecules or the degrada- 
tion of the RNA duplex is slow and ineffective, the 
recovery of target expression may be fast even when 
the sRNA regulator is fully recycled. 


Simulation results versus experimental results 


To test our model, we simulate the time evolutions 
of sRNA and its target mRNA, respectively, when 
they are expressed alone and together in the system. 
As shown in Fig. 8a, according to our simulation 
with reasonable parameters, more than half of 
mRNA molecules are inhibited after 1 h of expres- 
sion. For the fully recycling ratio 6=1, our simulation 
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Fig. 8. Time evolutions of sRNA and its target mRNA. 
First, SRNA was expressed over 10 h to arrive at steady 
state. Then at time 0, the sRNA gene was shut off, and 
mRNA was expressed for 65 min. Parameters are set as: 
as=4nM min‘, a =7 nM min‘, and k=0.001 (nM miny '. 
Other parameters are set as in Fig. 2. (a) Time evolution 
of sRNA and mRNA with ô=1. Superscript “+” or “—” 
indicates the presence or the absence of relevant RNA, 
respectively. (b) Effects of 5 on time evolutions. Empty 
symbols represent sRNA, and filled symbols represent 
mRNA. 


result shown in Fig. 8a is consistent with experimen- 
tal results." 

We focus on the time point where the number of 
mRNA and the number of sRNA are equal to each 
other in a cell. After that point, mRNA continues to 
accumulate, while sRNA is still actively degraded 
(see Fig. 8b, gray rhombus). Consistent with the 
previous analysis, the more effective is the recycling, 
the more time it takes to reach the time point to start 
gene expression. Under the same condition, it takes 
twice as long for target mRNA to recover from the 
state of suppression with null recycling (6=0) than 
with full recycling (6=1). 


Discussion 


MicM was first proven to be an unusual antisense 
sRNA when it did not destabilize even when target 
ybfM mRNA was overexpressed. However, for 
another target gene that belongs to a chb operon 
and possesses complementarity to the same region 
of MicM, annealing of MicM causes coupled 
degradation. For the most identified antisense 
sRNA-mediated gene silencing in bacteria, the Sm- 
like protein Hfq binds efficiently to both types of 
RNAs to stabilize them _'7-'? and to accelerate 
annealing between them.'”'® Here we start with a 
basic model of MicM-ybfM mRNA interaction and 
study in detail the noise properties of target 
expression by simulation. Our calculation shows 
that the inhibition effect of sRNA will be weak even 
with a high recycling rate (target mRNA will escape 
easily as k is reduced; data not shown). It implies 
that the inhibitory efficiency of catalytic sRNA 
depends more on the efficiency of sRNA—mRNA 
duplex formation than on the recycling ratio. On the 
other hand, it has been shown in Fig. 6b that a 
weaker binding between sRNA and mRNA usually 
causes more noise on target expression in the 
silencing region. All these findings may contribute 
to a better understanding of physiological signifi- 
cance and further design of the small RNA gene 
regulatory network. 

It is interesting to note that all curves of mRNA in 
Fig. 2 can be considered equivalent if we scale a, to 
#5 and if we scale \ to 74 In other words, the 
recycling ratio linearly affects the steady-state 
behavior of mRNA, but its influence on noise is 
distinctive; both the peak width and the peak value 
of the noise curve are significantly affected by the 
recycling ratio in a somewhat complicated manner, 
which will be considered in further studies. 

It is known that sRNA is prevalent in bacteria, 
and it acts by various mechanisms besides base- 
pairing with other RNAs. For example, a few of the 
sRNAs are part of ribonucleoproteins that carry out 
specific housekeeping functions. Some others act by 
mimicking the structures of other nucleic acids to 
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sequester their cognate proteins.* All these regula- 
tory modes are supposed to make sense in the 
selective evolution of life. Intricate regulatory infor- 
mation can be deciphered regardless of whether the 
regulator is a protein or an sRNA, the level of 
regulation, and whether the regulation is enzyme-like 
or not. In this sense, different noise performances 
may ensure life activity, diversity, and complexity. 


Materials and Methods 


The analysis was carried out using a simplified rate 
equation model and Monte Carlo (Gillespie algorithm”®) 
simulations. All simulation parameters were obtained 
from a reasonable estimation of experiments. 1> 

The Gillespie algorithm, developed by Daniel T. 
Gillespie in the 1970s, is well known as an effective SSA. 
By calculating the reaction probability density function 
P(t,) using the current reaction rate and reactant 
species combination number, we can determine the 
probability of when the next reaction will occur (T) and 
which reaction will occur (u). Then we renew the 
molecular number of relevant species. The Gillespie 
direct method needs two uniform random numbers to 
generate the random values r and pu that obey the 
reaction probability density function, so it is also a type 
of Monte Carlo algorithm. In this way, we obtain 
statistically correct trajectories of the chemical system's 
time evolution, which is rigorously equivalent to that 
produced by the grand probability function P(X1,...,Xn;t) 
as the solution of the master equation. 

It is quite difficult to obtain an integrated analytic 
solution for even a simple master equation. The most 
important advantage of SSA is that it is an exact algorithm 
that takes full account of fluctuations and is easy to code in 
a computer program. It is also effective for calculating 
other relevant physical variables such as the Fano factor 
and response time. 

Based on the ergodic hypothesis of equilibrium statis- 
tical physics, the ensemble average coincides with the time 
average (see Eq. (6)). To obtain the mean and variance of 
the molecules of a biochemical system, we only apply the 
frequency statistics (for a long-enough time) of the 
stochastic state sequence produced by SSA in steady 
state, with state waiting time as weight: 


T 
(%) = APC) = Jims f, AXO (6 
where X represents the value of any reaction species here, 
and the function f(X) may represent its mean value or 
variance. 

We use a so-called MFPT to describe the response 
characteristic of the biochemical system. MFPT is 
defined as:** 


a Thx; 
(Ta) = | t( - dP,), whereP; = XO P(Xi,..., Xyit) 
0 X;=0 
(7) 


Here, X; is the variable of interest, and Thy, is its 
threshold value, while (Tẹ represents its MFPT crossing 


the threshold. It theoretically represents the time it takes, 
on average, to reach a boundary the first time. The 
threshold values here mean certain boundary conditions. 
However, by simulation, we only use the stochastic state 
sequence produced by SSA and set a threshold value for 
the chemical reactant that we want to analyze and then 
calculate the time when it crosses the threshold the first 
time. Finally, we obtain the MFPT by averaging the 
results that are repeatedly derived. 
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